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Abstract 

In transport phenomena, perturbation waves are a result of interaction of molecules 
in gases and liquids, charged particles (ions, electrons) in plasma, conduction elec- 
trons and phonons in solid bodies. General statistical theory of the perturbation 
waves is developed and its corollaries are studied. On this basis is proved universal- 
ity of introduced earlier local time concept, which leads to a formulation of kinetic, 
conservation and governing equations for macroscopic transport phenomena with 
finite speed of the perturbations propagation in gases, liquids, solids and plasma. 

Speed of thermal perturbations propagation in phonon and Fermi electron gases 
and plasma, and also speed of thermal, momentum and mass perturbations propa- 
gation in ideal gas are theoretically determined. 

It is shown that published experimental results for femtosecond laser heating of 
thin gold films and results of power modulation experiments in JET tokamak agree 
with the developed theory. 

Keywords: Transport phenomena, perturbation waves, kinetic theory. 

1 Introduction 

It is well known that if a heat conducting system is initially in thermal equilibrium, and a 
local thermal perturbation is introduced, then according to Fourier equation temperature 
everywhere in the system instantaneously becomes perturbed. It means that according 
to Fourier equation, speed of thermal perturbations propagation is infinite. All other 
classical transport equations for viscous flow and diffusion also predict infinite speed 
of perturbations propagation. In reality, perturbation waves have finite speed, because 
perturbations propagation is a result of interaction of molecules in gases and liquids, 
charged particles (ions, electrons) in plasma, conduction electrons and phonons in solid 
bodies. 



In my previous works PQ-[3], approximate statistical theory of the perturbations waves 
in transport processes and phenomenological thermodynamic theory of thermal pertur- 
bations propagation were developed. They led to a proof of the local time concept for 
systems with zero or small macroscopic velocity, which leads to a formulation of kinetic, 
conservation and governing equations for real transport processes with finite speed of 
the perturbations propagation. Correctness of the local time concept was confirmed by 
analysis of experimental data for electron temperature pulses propagation in magnetically 
confined hot plasma. 

In the present work, general statistical theory of the perturbation waves is developed 
and its corollaries are studied. From the theory follows that the local time concept is 
universal and is applicable to all transport processes in non-equilibrium systems in gases, 
liquids, solids and plasma. 

I theoretically determine speed of thermal perturbations propagation in phonon and 
Fermi electron gases, and speed of thermal, momentum and mass perturbations propa- 
gation in ideal gas. For plasma in strong external magnetic fields, small values of the 
perturbations propagation speed in direction perpendicular to the field are theoretically 
predicted. I also study influence of macroscopic velocity on speed of the perturbations 
propagation. 

I show that published experimental results for femtosecond laser heating of thin gold 
films are in good agreement with the theory of the perturbation waves and the local 
time concept. I also find that experimental and calculated theoretical values of speed of 
thermal perturbations propagation in electron Fermi g close. 

I prove that published results of power modulation experiments with magnetically con- 
fined hot plasma in JET tokamak confirm correctness of the theory of the perturbation 
waves, the local time concept, and theoretical prediction of small values of the perturba- 
tions propagation speed in direction perpendicular to the field. These experiments also 
support theoretical prediction of equal values of thermal perturbations propagation speed 
in electron and ion plasma components. 

2 Statistical formulation of the problem 

I analyse a general case of non-equilibrium systems and use term particle for molecules 
in neutral gases and liquids, free ions and electrons in plasma, phonons and conduction 
electrons in solids. It is assumed that a non-steady, n on- equilibrium state field of thermo- 
dynamic and hydrodynamic parameters, for instance temperature, pressure, concentration 
and macroscopic velocity, exists initially in a region. Then some part of the initial field 
is perturbed. As a result, interacting perturbed and non-perturbed subregions appear 
in the region. Further analysis applies the same statistical method that was developed 
in my previous works PQ, [2] where I analysed a particular case when initial field is in 
equilibrium or steady non- equilibrium state. 
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Let E be an arbitrary point of the surface Sp separating the perturbed subregion 
from the non-perturbed one. At that point a local cartesian coordinate system X, Y, Z 
is introduced, where axis X is normal to the surface Sp at the point E, and axes Y 
and Z create a plane which is tangential to the surface Sp. The coordinate system 
X, Y, Z does not move with respect to the non-perturbed subregion. The axis X has 
direction from the perturbed subregion to the non-perturbed one. For the perturbed 
and the non-perturbed subregions distribution functions for selected kind of particles are 
designated as /(r, v,r) and /o(r, v, r), respectively, where r(X,Y,Z) and v(vx,vy,vz) 
denote the particle position vector and particle velocity vector, respectively; r is time. 
At the point E, particles with positive vx velocity component vx = > enter the 
non-perturbed subregion from the perturbed subregion, and particles with negative vx 
velocity component Vx — v%~ < enter the perturbed subregion from the non-perturbed 
one. If F is distribution function at the point E, there is F = /(r, v, r) for particles 
with vx = Vjp > 0, and F = / ( r > v ? r ) f° r particles with vx = < 0. Here and 
further superscript (+) is used for particles with positive vx = v% velocity component, 
and superscript (— ) - for particles with negative vx = velocity component. 

Let ips be a property of a particle crossing the surface Sp, ip(r,v,r) is the same 
property of a particle in the perturbed subregion, while y9 (r,v,r) is the same property 
in the non-perturbed subregion. Thus, <pg = (p(r, v,r) for particles with vx = > 0, 
and (ps — V^oty) when vx = < 0. 

The distribution function F satisfies general kinetic Boltzmann equation applied for 
selected kind of particles 

OF OF OF OF „ 

-^ + vx M + v Yw + v z -=D c -G M (1) 

where Dc is a collision term; Gm is a momentum term 

G _ dpx_dF_ + dp^dF_ ^dp^dF_ 
M dr dp x dr dp Y dr dp z 

and p{px,Py,Pz) is a particle momentum. We do not introduce any limiting assumption 
regarding structure of the collision and momentum terms in equation Q. 

As usually [1] , both sides of equation (jlj are multiplied by (p$ and after simple trans- 
formations it becomes 

d(<p s F) d(<p s F) dfoF) djysf) 

^^ +Vx ^t +Vy ^y^ +Vz ^z^~ 

= <p s {Do - G M ) + F(- W + vx^ + v Y — + v z —) (3) 
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All terms of expression ^ are integrated over velocities v, and integrals are trans- 
formed taking into account formulated earlier properties of functions F and (ps- It is sug- 
gested that for all values of velocity components v% , vy, fz derivatives , ^§p--, 
do not change the signs, and for all values of velocity components \ vy, vz derivatives 

9 %x°^ ' 9 %y oS> ; 9 %°z°' > a ^ so ^° no ^ change the signs. Therefore the integral mean value 
theorem [5] is applicable. After integration the following expressions are obtained: 

I T + I x + Iy + Iz = Icm + Is (4) 
d(<f S F) 



dr 



dvx dvy dvz 



dvv^ dvy dvz + [ ^lUL d v ( ) d Vy d V2 
or J or 



dr + dr [b) 
d(<f S F) 



Ix = J vx gx dv x dv Y dv z 



(+) f d(vf) , (+) , > . ~F) /" <%o./o) , (-) , , 
x ' dX x x J dX x 



T+) 9(n%W) -H- a^-W -)) 
~ + dX 



• ,r (6) 



J y = y ^y dv x dv Y dv z = 

= uy w y dv Y dv z + v Y ( ' j ^ dv x dv Y dv 2 



= ^ W ^7 + ^ j gf ( 7 ) 

Iz = J v z — — dv x dv Y dv z = 

— (+) f d ^f) A (+) A A , — (-) f d Mo) j (— ) j , 

= «z v ; y «^y «^ + ^ y dz «^y = 
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(8) 



Icm = J Vsi.Dc - G M ) dv x dv Y dv z (9) 

T f „ {dips . dips . dips . dip s \ , , , 

J \~dr~ Vx ~dX VY W + VZ TZ) dVx dVY dVz = 

= y ; [a; + ^ 9x + vy w + v ^ j dvY dvz+ 

= U \dr + Vx dX + VY W + Vz dz) +n [ Vx dX + VY W + VZ WJ (10) 

where n is mean number of particles per unit volume, overline denotes mean value of an 
appropriate parameter, and 

V& = Jv { X +)d ^ dv^ dvy dv Z j J M dvy dv Z = C X (11) 

vP=f V? dvtf dVy dv Z / J dvtf dVy dv Z = C 0X (12) 

W^ = J Vy 9 ^ dv^ dVy dv z j J ^ d V™ dVy dv Z = Cy (13) 

VP=J VY 9 ^ dv? dVy dv z j J dv? dVy dv Z = C y (14) 

vP=fvz?%p- dv^ dVy dv Z j J ^ dv^ dvy dv Z = C Z (15) 

— f d{p fo) (_) / fd(ip fo) , (-) , , , R v 

v z = J v z — — — dv x dv Y dv z / J — — — dv x dv Y dv z = c oz (16) 

denote respective mean integral values of v^\ Vjc\ Vy^\ Vy~\ v^~\ v z ^ By definition, 
c x > 0, c x < 0, and cy, c y, c z , c oz may have any signs or to be zero. 
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3 Evolution of the perturbed and non-perturbed sub- 
regions 



After substitution of formulas (|5j)-( 10 ) in Q we obtain expressions describing evolution 
of the perturbed and non-perturbed subregions 



Dcf) , D0 O 



+ 



Dt Dt 



A{t,X,Y,Z) 



(17) 



Dcf) 90 90 90 90 
D^ = d^ + Cx dX +CY W + Cz dZ 

D(f) 90o 90o 90o (900 

97 + Cox dx +C0Y W + Coz ~dZ 

A(t,X,Y,Z) = I cm + I s 



(18) 

(19) 
(20) 



where <f>(r,X,Y,Z) = n^+J, c x (r,X,F,Z) > 0, cy(t,X,Y,Z), c z (t,X,Y,Z) are 
parameters determined only by averaged properties of the perturbed subregion, while 
0o(r,X,Y,Z) 

only by averaged properties of the non-perturbed subregion. 



n ( Vo \ c ox(r,X,Y, Z) < 0, c y(r, X, F, Z), c z(r, X, Y, Z) are defined 



Macroscopic formulae (17)- (20) directly follow from the general kinetic Boltzmann 
equation (lj). They include total derivatives ^ and describing evolution of the 
boundary surface Sp separating the perturbed and non-perturbed subregions. As the 



perturbed subregion is in non-steady non-equilibrium state, therefore in formula (18) 
both local 1^ 7^ and convective cx §x + °y §y + c z ff 7^ derivatives are present. It 
means that total derivative ^ ^ describes propagation of the perturbed subregion into 
the non-perturbed one with finite total speed c(cx, cy, cz) and normal speed cx > 0. To 



prove it we multiply both sides of expression (18) by St which is a small time interval, 
and obtain following formula 



A0 = 0(r + St, X + c x St, Y + c Y St, Z + c z St) - 0(r, X, Y, Z) 



Dt 



St 



d~T +Cx dX +CY dY +Cz dz) 6T 



(21) 



Formula (21) shows that if X,Y,Z are coordinates of a point, belonging to the per- 



turbed subregion and located at the initial time moment r on the surface Sp, which 
separates the perturbed subregion from the non-perturbed one, at the next time moment 
t + St the field in a point with coordinates X + cx St, Y + cy St, Z + c z St becomes 
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perturbed, while at initial time moment r it was not perturbed. It means that the bound- 
ary of the perturbed subregion gradually propagates in the non-perturbed subregion with 
finite normal mean speed cx{t, X, Y, Z). Therefore cx{t, X, Y, Z) is local mean speed of 
the perturbations propagation. 

In the same way we can prove that if the non-perturbed subregion is in non-steady non- 
equilibrium state, the total derivative 7^ describes propagation of the non-perturbed 
subregion into the perturbed one with finite total speed c (cox, cqy, coz) and normal speed 
Cox < 0. Therefore in this case, at initial moment of time, when the perturbation was 
introduced, the boundary surface Sp instantaneously disappears, two subregions - perturbed 
and non-perturbed - join, and the perturbations propagation is principally non-observable. 

If the non-perturbed subregion is in steady non-equilibrium state, local derivative 
disappears in formula (19), and only convective derivative cqx §^ + C ° Y Ity + 

derivative does not describe propagation of the 



9<j>o 
dr 



Now the tota 



c oz H ^ is left. 

non-perturbed steady state subregion into the perturbed one, because such propagation 
must be accompanied by change of respective properties in time, and propagation of 
steady state field is physically impossible. In a particular case of equilibrium in the 
non-perturbed subregion there is = 0. 

Therefore we proved that in all non- equilibrium processes the perturbations propaga- 
tion is observable and physically meaningful only when the perturbed subregion interacts 
with steady state non-perturbed subregion, and only propagation of the non-steady state 
perturbed subregion into the steady state non- equilibrium or equilibrium non-perturbed sub- 
region is physically possible. It means, that in all cases of gas, liquid, solid and plasma 
non- equilibrium systems, the non- equilibrium perturbed subregion of any three-dimensional 
shape gradually replaces the non-perturbed steady state subregion with local normal speed 



cx = > which is speed of the perturbations propagation defined by expression (11). 



By definition, speed of the perturbations propagation cx > is positive and finite. 
The limitations introduced by the integral mean value theorem [5] ensure that for any 
arbitrary distribution function / and particle property ip, expression (11) determines 



positive, finite and physically sound values of speed of the perturbations propagation. If 



these limitations do not hold, the expression ( 11 ) does not define positive and finite values 



of cx, and does not serve more as statistical definition of the perturbations propagation 
speed. In such case, other methods of cx determination can be applied, for instance, 
numerical statistical modelling of the perturbed subregion propagation or experimental 
study of the perturbations propagation. When the integral mean value theorem does not 
hold, approximate values of speed of the perturbations propagation can be calculated if 
instead of derivative ^feP to use its absolute value 



; (+) 
J x 



; (+) 
J x 



dX 



dv^ dvy dv 



Y 



1 gfefl I 
ox I 

d(<pf) 



dX 



dvjp dv Y dvz 



(22) 



When does not change its sign during integration, formula (|22|) is identical to 



(11) and defines theoretically exact value of cx, in the opposite case it defines approx- 
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imate value of speed of the perturbations propagation, whose accuracy is checked later 
by comparison with experimental data for thermal perturbations propagation in Fermi 



electron gas. Thus, expression (22) is universal, and we will always use it for theoretical 



determination of speed of the perturbations propagation. 

4 Perturbation traveltime and generalization of the 
local time concept 

As in my previous works PQ-[H], I assume that perturbation traveltime Tp(x,y,z) defines 
a time moment when the perturbation reached a given point M(x, y, z), where x, y, z is 
an arbitrary cartesian coordinate system not moving with respect to the non-perturbed 
steady state subregion. Therefore the boundary surface Sp, separating the perturbed from 
the non-perturbed steady state subregion, is a surface of constant perturbation traveltime 
Tp(x, y, z) = const. From definition of perturbation traveltime follows 

|Vr P | = — r (23) 

c x {Tp,x,y,z) 



or 



dr P \ 2 ( dr P \ 2 (dr A 



2 



In the equation (24) instead of speed cx(Tp,x,y, z) can be used slowness 



s x (r P , x, y, z) = — - — 1 r (25) 

c x {Tp,x,y,z) 



The non-linear equation (24) is similar to the classical eikonal equation [6], [7] and 
determines traveltime in the case of any macroscopic velocity. According to results ob- 
tained in previous sections, only its solutions corresponding to primary perturbation waves 
propagating in the non-perturbed steady state subregion, have physical meaning. 



From equation (24) follows that 

c x (T P ,x,y,z) = c x (r P (x,y,z),x,y,z) = c x (x,y,z) (26) 

s x (r P ,x,y,z) = s x (r P (x,y,z),x,y,z) = s x (x,y,z) (27) 
are functions only of coordinates. 
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In a general case, the non-perturbed steady state subregion can be perturbed in differ- 
ent places and at different time moments. Accordingly solution of eikonal equation (24) 
must satisfy several initial conditions of the type 



T 0P (x,y, z) = T 



(28) 



where T p(x,y, z) and r denote initial value of traveltime and time moment when the 
perturbation was introduced, respectively. 

So if the non-perturbed initial steady state subregion was perturbed in different places 
and at different time moments, perturbed non-steady state subregion is permanently divided 



on different sub-subregions corresponding to different initial conditions of type (28). 

For any point M(x, y, z) of the region, local time is determined as time counted from 
a moment when the perturbation reached the point p]-[3] 



$(x,y,z,T) = t - T P (x,y,z) 
For this point at a global time moment r, three characteristic cases exist: 



(29) 



1. i? = t — Tp < 0, the perturbation has not reached the point. Thus, the point belongs 
to the non-perturbed steady state subregion. 

2. d = t — Tp = 0, the perturbation has reached the point M. Now M is located on 
the border surface Sp, separating the perturbed and the non-perturbed subregions. 

3. $ = t — Tp > 0, the point M is inside the perturbed subregion, where transport pro- 
cesses occur. All further considerations are formulated for the perturbed subregion 
where $ > 0. 



For a case when the perturbed subregion interacts with steady state non-perturbed 



subregion, fundamental equation (17) has a form 



9<f> , d<j> d(p dcp 

^ + Cx dX +CY dY +Cz dZ ==MT > X > Y > Z) 



(30) 



where A = A — (c x + c y ^ + c z t^f)- In a case of non-zero macroscopic velocity 



there is cy ^ 0, cz ^ in expression (30) 



As 



Further in formula (30 ) we can replace independent variables r, X, Y, Z by X, Y, Z. 



d<p dcj) d<p 
~dr~ = <9tf ' dX 



OX 



« c x ft?' dY 



d(f) 


d(f) 


d(f) 


T ~ dY 


*' dZ 


T ~ dZ 
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equation (30) becomes 



cx 



dX 



+ Cy 



dY 



oZ 



c ■ V0 = A c 



(31) 



where symbols | r , 
respectively. 

Finally we get 



mean that independent variables are r, X, Y, Z or X, Y, Z, 



dX r 



An 



(32) 



where coordinate X c is measured along axis parallel to vector c(cx, cy, cz), and c 



Cx + Cv + c 



Y 



Left side of the general fundamental equation (30) includes first order time and space 



derivatives and explicitly depends on velocity components cx,cy,cz- Expression (32) is 



fully equivalent to equation (30), but its left side does not depend explicitly on cx, cy, cz, 
and speed of the perturbations propagation cx affects only local time Because of this, 



transformed equation (32) is identical for finite cx, when Tp > 0, and for the classical 
case of infinite cx, when always Tp = 0. The only difference is: for finite cx there is 
< d < r, while for the classical case always d = t. Therefore when independent 
variables X, Y, Z are used in all transport equations, the same equations describe 
the classical case of infinite speed of the perturbations propagation when i? = r, and the 
case of finite speed of the perturbations propagation when < i? < r. Thus, we proved 
that the local time concept, formulated in my previous works [l]-[3] for systems with zero 
macroscopic velocities cy = 0, cz = 0, is universal and holds for a general case of non-zero 
macroscopic velocities when cy ^ 0, Cz ^ 0. It means that: 

For all non- equilibrium systems in gases, liquids, solids and plasma, any macroscopic 
velocities and all transport processes, kinetic, conservation and governing equations in a 
case of finite speed of the perturbations propagation and the respective classical equations 
with infinite speed of the perturbations propagation (for instance, Fourier, Fick, Euler, 
Navier-Stokes and magnetohydrodynamic equations) are identical , if independent variables 
d = t — Tp, x, y, z are used. The same equations become classical when Tp = 0, •& — r 
that corresponds to cx — > oo, and they describe the finite cx case when t p > 0, < d < t. 
Because of this, all equations written in such a form are called modified equations.. 

For classical transport processes such as heat conduction, diffusion and viscous flow, 
modified equations and some analytical solutions of modified governing equations are 
presented in papers pQ- [3]. 

Immediate corollaries of the local time concept are [XJ - [3J : 

1. Modified governing equations for all transport processes with finite speed of the 
perturbations propagation are of the same parabolic type, as classical equations, 
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thus, fully excluding any possibility of reversing the transport process. These equa- 
tions are consistent with the Second Law of thermodynamics because their solutions 
satisfy the Maximum Principle. According to the Maximum Principle for parabolic 
partial differential equations without a source term, maximum value of any scalar 
characteristic parameter, e.g. temperature, is reached either on boundary of the 
domain or at the initial moment of time. 

Numerous earlier attempts of many authors to introduce transport equations with 
finite speed of the perturbations propagation were not successful because they tried 
to formulate hyperbolic governing transport equations which are not consistent with 
thermodynamics, and therefore are non-physical. For instance, some their solutions 
describe flow of heat from cold to hot bodies in contradiction with the Second Law 
of thermodynamics. Full enough list of these works and their analysis I give in 
papers PQ-[3]- 

2. Let 0(x,y, z,t) be a solution of a classical transport governing equation, satisfying 
certain boundary conditions and an initial condition (3(x, y, z, 0) = (3 (x, y, z), where 
f3 is a scalar characteristic parameter. According to the local time concept, a solution 
of the modified governing transport equation is the same function but with local 
time as an argument /3(x, y, z, $). 

Solutions of classical parabolic type governing transport equations predict that any 
local perturbation introduced at initial moment of time Tj = 0, instantaneously 
affects all space domain. The modified parabolic governing equations use local 
time $ as an independent variable. In this case, initial perturbation introduced at 
initial local time moment t?j = 0, affects the space domain at the same local time 
value i9j = r im — Tp(x,y,z) = 0. Therefore perturbation arrives to an arbitrary 
point M(x,y,z) at global time moment T im = Tp(x,y,z), i.e. with global time 
delay. So the modified parabolic governing equations predict finite speed of the 
perturbations propagation, and the perturbations have wave behaviour. As a result 
of transport processes irreversibility, these waves cannot reflect and interfere because 



only solutions of eikonal equation (24) corresponding to primary perturbation waves 
propagating in the non-perturbed subregion, have physical meaning. 

On the boundary of the perturbed subregion it is $ = 0, and therefore on the bound- 
ary is (3(x, y, z, i?) = /3q(x, y, z). Therefore the propagating wave of any characteristic 
parameter does not have a sharp front. 

From formulae derived in previous paragraph follows that solutions of classical 
transport governing equations are accurate enough when the transport process is 
slow and/or perturbation traveltime is small enough. To prove this, we choose 
an arbitrary point with coordinates xa, yA, Z A an d perturbation traveltime value 
tpa(xa,Va,za), where there is 0(x A , Va, za, $ = r - t pa ) = (3 (x,y,z) at initial 
moment of time r < tpa- Then becomes r > tpa, and {3(xa, yA, za, $ = t — tpa) ^ 
Po(x,y,z). Obviously, the difference between solutions of classical and modified 
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governing transport equations for the perturbed subregion 



| = \/3(x A , Ua, za, t) - f3(x A , Va, z A ,d = t - t pa ) 



dp 



Or 



TPA 



(33) 



is small for small values of perturbation traveltime Tp A and/or slow processes, when 



1 1^ I is small. 

1 Or 1 



5 Modified heat conduction equations 

As an illustration, the modified Fourier kinetic, conservation and governing equations 
for heat conduction with finite speed of heat propagation in isotropic continuum are 
reproduced from [I]-[3] 



-K VT 



S(x,y,z,$), 



(34) 
(35) 



1 dT 



V 2 T 



(36) 



where q, T, K, C denote heat flux, temperature, heat conductivity, and specific heat 
capacity per unit volume of the continuum, respectively; S(x, y, z, ■#) is a heat source 
term; k = K/C is thermal diffusivity. It is assumed that K = const. For infinite speed 



of the perturbations propagation, when $ = r, equations (34)-(36) become classical. 



For regular independent variables r, x, y, z and finite speed of the perturbations 



propagation, equations (34)-(36) are transformed into 



dT 

q = —K ( VT + — ■ Vr P 



(37) 



dT da 
C -^ + V-q+^-VT P = S(x,y,z,$), 



(38) 



Equations (37), (38) and (39) are identical to (34), (35) and (36), respectively, and are 



quite different from classical Fourier equations predicting infinite speed of propagation. 
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Figure V. Dimensionless temperature distribution from instantaneous plane heat source 
(1 - according to modified Fourier equation, 2 - according to classical Fourier equation). 



Only when cx — > oo, and therefore Tp — > 0, they become classical. We see that using 
regular variables r, x, y, z instead of x, y, z, makes equations describing transport 
process with finite speed of the perturbations propagations, much more complicated and 
more difficult to solve. 

In one dimensional case and for zero heat source term modified Fourier governing 



equation (36) becomes 



dT 



d 2 T 



dx" 2 



(40) 



x 



We assume initial condition T = for d = and instantaneous source Q at plane 
= and $ = 0, while Q C is quantity of heat liberated per unit area of the plane. For 



this case solution of equation (40) is 



T 



Q 



exp 



x 

4^9 



(41) 



Expression (41) is obtained by applying the local time concept to a respective solution of 



regular Fourier equation presented in [8]. 
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As d = r — Tp = r—\x\/cx, t = r d K/c x , \x\ = x d n/cx, where r d , x d are dimensionless 



time and dimensionless absolute value of space coordinate, respectively, formula (41) can 
be transformed into 



T r d 

7p = \ ex P 

T V T d - %d 



x: 



4(r d - x d ) 



(42) 



where Tq is temperature at x — 0. 



Curve 1 in Fig. 1 is temperature distribution according to expression (42) for r d = 1, 
and curve 2 is temperature distribution according to classical Fourier equation for the same 
value of r d = 1. We see that modified Fourier equation really describes heat conduction 
with finite speed of the perturbations propagation. 



6 Theoretical determination of speed of the pertur- 
bations propagation in transport processes 

In this section, we theoretically determine speed of the perturbations propagation cx in 
ideal gas, solid and plasma. It is assumed that ideal gas, Fermi electron gas in solid, 
and electrons and ions in plasma are in very close to equilibrium state. As a result, 
equilibrium distribution functions can be applied for calculations and analysis. It can 
be proven that in this case there is \d((pf)/dX\ = A C (X,Y, Z) B(T,vx,vy,vz) where 



T is temperature. Therefore formula (22) for calculating speed of the perturbations 



propagation is transformed into an expression not containing space derivatives 

M +) Bdv^ dVydVz 

cx = rrj (43) 

J B dv x dvy dvz 



system is zero = 0, expression (43) becomes simpler 



For a case when B(T, t>), where v = \ v\ + v Y + Vg, and macroscopic velocity of the 



/ v 3 B(v) dv , AA . 

Cx = O f 2 Of \ a 44 
2 J v 1 B[v) dv 

6.1 Speed of the perturbations propagation in a moving media. 



Expression (|22j) defines speed of the perturbations propagation in a general case. For a 
moving media, particles velocity v is a sum 

v = v r + v h (45) 
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where Vh( v hX,VhY,Vhz) is macroscopic velocity, in a case of fluid it is hydrodynamic 
velocity, and v r (v r x, v r y, v r z) is particles velocity in non-moving media when = 0. 



From expressions ( 22 ) and ( 45 ) follows that 



c x = VhX + 



(46) 



poo roo p 

h= y y 

J v rX =-VhX Jv rY =-oo Jv rZ =-oo 



V r X 



OX 



dv r x dvrY dv 



hZ 



v rX = -VhX Jv rY =-oo Jv rZ 



dX 



dv r x dv r y dv 



rZ 



(47) 



(48) 



where it is assumed that the distribution function / depends only on v r (v r x, v r y, v r z) 
velocity components. 



Expressions (46)- (48) show that speed of the perturbations propagation increases when 



macroscopic velocity component v^x increases, and becomes zero only in a non-realistic 
case of Vhx — > — oo. The function cxivhx) has two asymptotes: at Vhx > a line cx = Vhx 
is an asymptote, and at Vhx < a line cx = is the second asymptote. At small Vhx 
there is 



cx ~ cqx + k c v hX 



(49) 



where c x is speed of the perturbations propagations at Vhx = 0, and 1 > k c > 0. Formula 



(49) shows that even for small values of macroscopic velocity, speed of the perturbations 



propagation in a moving media is not an additive function of the speed of the perturbations 
propagation in non-moving media and macroscopic velocity. 

All these results regarding function Cxiyhx) are obtained by using only two assump- 
tions: (a) the improper integrals I a , /& converge, and (b) a function 



Ic(Vrx) 



oo poo 



— CO J —CO 



d(<pf) 



dX 



dv r Y dv. 



rZ 



(50) 



does not depend on sign of particle velocity component v r x- 



6.2 Isothermal viscous flow of ideal gas 

Speed of the perturbations propagation in isothermal ideal gas viscous flow is calculated 
using formulae (46)- (48) and Maxwell distribution function from [1] 



/ 



n 



m 
2M 



cxp 



mvj\ 
2kf) 



(51) 
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Figure 2: Dimensionless speed of the perturbations propagation Cx as a function of di- 
mensionless hydrodynamic velocity normal component Vhx for ideal gas isothermal viscous 
flow (1 - Cx-, 2 - asymptote). 

where n, m, k denote number of particles per unit volume, particle mass and the Boltz- 
mann constant, respectively; v% = v^ x + v%y + v rZ- 

For viscous flow, the property of interest ip is a vector of tangential momentum of a 
particle 



where j, k are unit vectors parallel axes Y and Z, respectively As Vh Y , Vhz and gas 
density p = n m are functions of space coordinates X, Y, Z and time r, we get 



It is convenient to normalize all velocities using most probable particle speed vq = 
(2/cT/m) - 5 from [4j as a reference: C x = c x /v , V hX = v hX /v , V rX = v rX /v , V rY = 
v rY /v , V r z = v rZ /v , V r = v r /v . 



ip = m {v hY + v rY )j + (v hZ + v rZ )k 



(52) 




(53) 
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From expressions (43), (46)-(48) and (53) follows that dimensionless speed of the 
perturbations propagation in isothermal ideal gas viscous flow is 



Cx = V hx + I f = V hx+ -L eX f~f\ (54) 
h V 71 " erfc{-V hX ) 



where 



poo roo poo 

h= / / Vrx exp(-V r 2 ) dV rX dV rY dV rZ = 

JV rX =-V hx JV rY =-oo JV rZ =-oo 

= |exp(-y ft 2 x ) (55) 



rOO roo roo 

/ / / exp(-V r 2 ) dVrx dV rY dV rZ 

JVrx=—Vh.x JV rY =-oo JV r z=-oo 



■K\/7r 



erfc(-V hX ) (56) 



Formula (54) and Fig. 2 confirm all previous general remarks regarding dependence 
of speed of the perturbations propagation from macroscopic velocity. For Vhx = speed 
of the perturbations propagation in isothermal viscous flow is simply 



1 2kT 2kT . . 

c x = ^\ = 0.56419 \ 57 



6.3 Conduction of heat in ideal gas 

Speed of the perturbations propagation for a case of heat conduction in non-moving ideal 
gas is determined for an arbitrary temperature distribution of T(X, Y, Z, r) and constant 
pressure conditions p = const. As n = p/kT, the distribution function (plj) becomes 



/=p(£)W«p(-^) (58) 

The property of interest for conduction of heat is kinetic energy of a particle tp = 
mVr/2. After simple calculations and by using introduced in the previous subsection 
dimensionless quantities we find that 

d Ml = A c (X,Y,Z) B(V r ) 
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where 



8T 

A C (X,Y,Z) = 2.5(m/2nk) 1 - 5 T- 2 - 5 — , B(V r ) = V 2 (0.4 V 2 - 1) ex V (-V 2 ) 



It means that the integral mean value theorem does not hold, and formula (44) is approx- 
imate. Numerical integration gives 

If 

C x = ^y- « 0.81755 (59) 

2 lg 

POO 

I f = V r 5 |0.4 V 2 - 1| exp{-V 2 ) dV r « 5.30405 (60) 



roo 

I g = V; 4 |0.4 V 2 - 1| exp(-1/ 2 ) dV r « 3.24387 (61) 

w 

and 



2kT 

c x ~ 0.81755 W (62) 

V m 

6.4 Isothermal self-diffusion in ideal gas 

For determination of speed of the perturbations propagation for isothermal self-diffusion 
in non-moving ideal gas we assume an arbitrary space distribution of labeled particles 
quantity per unit volume n(X, Y, Z) at T — const. As for labeled particles Maxwell 
distribution function holds, and property of interest is the particle velocity <p = v r , we 
obtain 

d((pf) dn f m \l / mv 2 \ 
dX =Vr dX\2nkT / 



Therefore the mean integral value theorem and relation (44) are applicable, and di 



mensionless speed of the perturbations propagation is easily calculated 

C x = ^r = 0.66467 (64) 



where 



roc 1 /5\ 

I h = J Q V* exp(-V 2 ) dV r = - r|-J (65) 



18 



poo I 

I, = / V r 3 exp(-V r 2 ) dV r = - T(2) 



(66) 



Finally we get 



2kT , s 

c x = 0.66467 \ 67 

V m 



Integrals (65) and 66) are calculated using method described in [9]. 

We see that speed of the perturbations propagation in ideal gas is different depending 
on type of transport process: at the same gas temperature, it is the highest for energy 
transport (thermal conductivity) at constant pressure, it is the lowest for momentum 
transport in isothermal viscous flow at Vhx — 0, and it has intermediate value for mass 
transfer in isothermal self-diffusion. 



6.5 Conduction of heat in solid 

In dielectric crystals thermal perturbations propagate due to phonons, and because of 
this speed of the perturbations propagation is equal to speed of sound in the crystal if we 
take into account only acoustic phonons [9]. In metals thermal perturbations propagate 
in Fermi electron gas much faster then in lattice, because propagation speed in Fermi 
electron gas has an order of magnitude of Fermi velocity. Naturally, propagating thermal 
perturbations in electron gas cause induced local thermal perturbations in lattice. As 
a result, speed of the thermal perturbations propagation in metals is equal to speed of 
thermal perturbations propagation in Fermi electron gas, and the perturbations traveltime 
space distribution is the same in electron gas and lattice. 

Fermi-Dirac distribution for electron gas at temperatures much lower than Fermi tem- 
perature is 

, 2m 3 p 1 . , 

/ = -nr — 1 (68) 



h 3 

1 exp 



r F \V 2 



+ 1 



where V = v/vp, rp = T/Tp, v, vf, m e , h, Tp denote dimensionless velocity, dimen- 
sionless temperature, electron velocity, Fermi velocity, mass of electron, Plank constant 
and Fermi temperature, respectively. 



Formula (68) is a modified version of classical expression presented in many sources, 
including for instance classical Born's book |10j . where is used Fermi-Dirac distribution 
function f e such that dn e = f e y/e de (n e and e are number of electrons per unit volume 
and electron kinetic energy, respectively), while we use distribution function / satisfying 
a condition dn e = f dv x dv y dv z . 
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Figure 3: Dimensionless speed of the perturbations propagation in electron gas Cx and 
gold temperature as functions of dimensionless temperature r F . 

For heat conduction in electron gas, parameter of interest is <p = m e v 2 /2, so 
d(<pf) 



OX 
where 



A C (X,Y,Z) B(V) 



A C (X,Y,Z) 



m\v\T F dT 
h 3 T 2 dX 



, B(V) 



V 2 (V 2 - l)exp 



T F \V 2 -1) 



expfr^ 1 ^ 2 -!)) +1 



It means that the integral mean value theorem is not applicable, and formula (44) is 
approximate 



C 



cx 



v F 2 J, 

h = V 5 \V 2 - 1| exp [r F \V 2 - 1) 



exp [r F \V 2 -l) +1 



dV 



(69) 
(70) 
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Il 



V A \V 2 - 1| exp (r F \V 2 - 1) 



exp r F 



^(V 2 -!)) +1 



dV 



(71) 



Fig. 3 presents values of Cx obtained by numerical integration of expressions (70) 
and (71). It shows that in wide range of gold temperatures from 300 K to 1100 K 



speed of thermal perturbations propagation, practically, does not change, and its average 
value is cx = 0.5003 Vp. As lim r . F _ s . Cx = 0.5, we can conclude that speed of thermal 
perturbations propagation in Fermi gas and metals does not vary significantly in wide 
range of temperatures - from very low to moderate high. 



6.6 Conduction of heat in plasma 

Speed of the perturbations propagation in plasma is different depending on whether exter- 
nal magnetic field is present or not. If the external magnetic field is absent and plasma is 
close to equilibrium, Maxwell distribution can be applied for our calculations, and speed 
of the thermal perturbations propagation in plasma electron component is determined by 



approximate formula (62) where must be introduced mass of electron m = m e . Mass of 



ions is much greater then mass of electrons 3> m e . However if we substitute m = rrii in 



formula (62) we do not receive correct value of speed of the thermal perturbations prop- 
agation in plasma ion component because thermal perturbations, propagating in plasma 
electron component, cause induced local thermal perturbations in the ion component. 
As a result, speed of the perturbations propagation in the ion component is equal to 
speed of the perturbations propagation in the electron component, and the perturbations 
traveltime space distributions are the same in both plasma components. 

When external magnetic field B exists in plasma, the most interesting is a case of the 
perturbations propagation in direction perpendicular to the field, because conduction of 
heat in this direction strongly affects efficiency of nuclear fusion devices such as tokamaks 
and stellarators with magnetic confinement of plasma. In strong enough magnetic field, 
gyrofrequencies of electrons and ions become much greater than their respective collisions 
frequencies. In this case, from qualitative considerations presented in |11J follows that 
common for plasma electron and ion components speed of the thermal perturbations 
propagation perpendicular to the field is 

cx ~ ~ (72) 

-DO 

So we can expect that in plasma with strong enough external magnetic field, speed of the 
thermal perturbations propagation perpendicular to the field is small. Experimental data 
presented in Ref. pQ supports this conclusion. 
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Figure 4: Thermal pulses traveltime vs. gold films thickness according to works [12] and 
[T5] presenting experimental results for femtosecond laser heating of thin gold films. 

7 Experimental verification of the perturbations prop- 
agation theory and the local time concept 

In this section, the perturbations propagation theory and the local time concept are 
tested using experimental data for conduction of heat in thin gold films and in hot fusion 
magnetically confined plasma. 

7.1 Femtosecond laser heating of thin gold films 

Experimental data presented in works [12] and [T3] makes possible to check validity of the 
local time concept for conduction of heat in metals, and also accuracy of our theoretical 
calculations relating to speed of the thermal perturbations propagation in metals. Authors 
of these works studied propagation of heat in thin gold films whose front surface was 
irradiated by very short femtosecond dye laser pulses with wavelength Ao = 630 nm, pulse 
duration was t\ p = 96 fs. They measured and registered the films back surface reflectivity 
history during and after irradiation. As the reflectivity change is caused by change of 
electron temperature, the cited experiments describe peculiarities of heat conduction at 
extremely short time intervals. Gold films thickness L varied from 500 to 3,000 A and 
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is significantly greater than gold optical skin depth do ~ 150 A. In experiments, focal 
spot diameter of laser radiation on the film front surface is d s ~ 20, 000 A>L, therefore 
quasi one-dimensional temperature T(x) and perturbation traveltime t p (x) distributions 
in a cylindrical control volume of diameter d c <C d s below the focal spot may be assumed, 
where < x < L is a space coordinate along the control volume. 

According to cited works [12] and [13] , change of films back surface reflectivity, which 
depends on electron temperature, starts with time delay after start of laser irradiation. 
It means that authors directly measured traveltime t p = of the thermal perturba- 
tion. Their experimental data shows that it approximately linearly depends on the film 
thickness L (Fig. 4). Using best fit (least square) approximation we calculated exper- 
imental dimensionless speed of the thermal perturbations propagation Cx = c x jvp = 
(vp dr p /dL)~ l = 0.6 ± 0.1 and film thickness Lq = 146 A corresponding to t p = 0. Thus, 
our theoretical value of Cx = 0.5 is in good agreement with experimental value of Cx- 
We see that L « do, it means that in the region < x < do there is r p = 0. This result 
is physically sound. 

Experimental data presented in [12J and [T3] makes possible to check the local time con- 
cept. For this purpose we will use a system of coupled one- dimensional modified Fourier 
governing equations describing electron gas and lattice thermal interaction according to 
the local time concept 

QT P d 2 T e 

Ke^r4- + 9e (Ti - T e ) + a (73) 



dd dx 

<?.-■?>> (74) 

where indices e, I denote electron gas and lattice, respectively; k(T) = K(T)/C(T) is 
thermal diffusivity; C(T), K(T) are heat capacity per unit volume and thermal con- 
ductivity, respectively; g = G/C(T); G is electron-phonon coupling constant; a = 
Ai($,x)/C e (T e ), Ai($,x) is a heat source term resulting from laser beam absorption by 
electron gas. 

For infinite speed of thermal perturbations propagation there is $ = r, and equations 



(73) and (74) become regular two temperature model equations from [H] with k\ = 0, 
obviously not applicable in analysed case because finite speed of thermal perturbations 
propagation was observed in experiments. 



The system of equations (73)-(74) is nonlinear. Its numerical solution cannot give 
clear answer regarding applicability of the local time concept, because the works [12] and 
[13] do not present absolute values of the films back surface reflectivity which is affected 
by electron temperature. The only available experimental data, besides traveltime, is 
back surface reflectivity in arbitrary units vs. time counted from zero-delay point of laser 
beam. From these data we can determine only time moment tq(L) when extremum of 
reflectivity dR/dr = is observed. Since R(T e ) is a monotonic function dR/dT e ^ 0, and 
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Figure 5: Applying the local time concept (diamonds) and classical concept (circles) to 
experimental data from [T2] and [13] for femtosecond laser heating of thin gold films. 



= r - r p (L), we obtain dR/dr = dR/dT e dT e /dr = dR/dT e dT e /dd = 0. Therefore 
extremum of reflectivity corresponds to maximum of temperature at local time moment 
$o(L) = Tq(L) — T p (L). So we can check the local time concept, if we will find analytical 
expression for $o(L) from the system of equations (73)- (74). 

In cited experiments the laser energy fluence was relatively small, approximately 1 
mJ/cm 2 . It results in moderate growth of the films temperature during laser heating and 
moderate changes of electron gas and lattice heat capacities. This allows to replace in 
equations (73)-(74) instantaneous values of C e , Ci, K e , ni by their average values without 
significant loss of accuracy. In this case, equations become linear with constant coeffi- 
cients, and we can derive their analytical solutions T e (x,$), Ti(x,$) corresponding to 
initial T e (x, 0) = Ti(x, 0) = T and boundary conditions dT/dx(0,$) = dT/dx(L,-d) = 
which we can assume for analysed case of laser heating. After introducing excess temper- 
atures t e — T e — T , ti = Ti — T and applying separation of variables we get 



t e = Y. X n{$) C0S[ —X 
n=l V L 



tl = J2 Y M C08[—3- 
n=l V L 



nn 



(75) 
(76) 
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In cited laser experiments with thin gold films there was C e <C Ci, K e 3> ki, g e ^> 
L/do > 3 and pulse duration rj p < 1 ps. For this conditions and an additional asymp- 
totic condition g e <C k\ the following asymptotic formula relating to electron temperature 



maximum is obtained from expressions (75)-(82) 



K, 



Ki 



'K 



A K z K + g e 
1 

ML) - ^oo 
1 

1 L2 



A K = IT 2 K e 



(83) 
(84) 

(85) 
(86) 



where "#oo < t~i v and 9j are parameters which depend on laser pulse duration and do not 
depend on the laser beam intensity. 



Formulae (83)-(86) are corollaries of the local time concept. Expressions (83)-(85) 



include known from experiment values of $q(L) and L, and four unknown parameters 
Ak, d e , ^oo, 9d- They are determined by least squares s = J2[Ki — ( y AKZK+d e )/9d} 2 = min 
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after introducing experimental values of -&q(L), L in (83)- (85). Then can be calculated 
n e = A k /tt 2 , T e = G/gtf, r K = n e / K es where G = (3.0 ± 0.5) ■ 10 16 W/m 3 K, 7 = 
66 J/m 3 K 2 [15j, and thermal diffusivity of gold is K es = 1.57 • 10~ 2 m 2 /s. 

To make a comparison, the same calculational procedure can be applied for a classical 
case of infinite speed of the perturbations propagation if instead of experimental local 
time values $o(L) experimental global time values r (L) are used. 

Fig. 5 presents experimental and best fit values Ki(Zk) for local time $o(L) and global 
time tq(L) coordinates of temperature maximums. Results of calculations using the lo- 
cal time concept are: s = 2.235 l/ps 2 , A K = 1.144 • 10 7 k 2 /ps, g e = 0.917 1/ps, 6 d = 
0.709, ^00 = 0.0715 ps, n e = 1.160- 10- 2 m 2 /s, T e = 496 K, r R = 0.74. Results of calcula- 
tions assuming infinite speed of the perturbations propagation are: s = 3.805 l/ps 2 , A K = 
3.486-10 6 A 2 /ps, g e = 1.730 l/ps, 6 d = 0.918, tf o = 0.0715 ps, K e = 3.532- 10~ 3 m 2 /s, T e = 
263 K, r K = 0.23. For a classical case was used the same value of $00 = 0.0715 ps as 
was determined by local time data because according to expressions (83)- (85) -#oo = 



lim L _>.o $o(L)- We could not calculate experimental errors because authors of cited works 
did not provide enough data. 

We see that experimental data from [12] and [13] for femtosecond laser heating of 
thin gold films support the local time concept and reject classical approach. In the 
last case, variance s is significantly greater, calculated from experimental data average 
electron temperature T e = 263 K is too low taking into account that all experimental 
equipment had ambient initial temperature, and determined from experiment electron 
gas thermal diffusivity is too small. Application of the local time concept gave value 
of average electron temperature T e = 496 K that is reasonably moderate and supports 



reasoning led to linearization of equations (73)-(74), while obtained from femtosecond 
laser experiments value of electron gas thermal diffusivity K e = 1.160 ■ 10~ 2 m 2 /s is close 
enough to n es = 1.57 • 10~ 2 m 2 /s. 

So we can conclude that all available experimental data for femtosecond laser heating 
of thin gold films presented in [12] and [13] support formulated earlier theoretical predic- 
tions relating to speed of thermal perturbations propagation in metals and the local time 
concept. 



7.2 Power modulation experiments with JET tokamak 

In my work [Tj the early versions of the perturbations propagation theory and the local 
time concept were verified for hot plasma by analyzis of experimental data for thermal 
pulses propagation in Alcator C-Mod tokamak. Now we will check validity of developed 
now theory and the local time concept for hot plasma in JET tokamak by analysing 
experimental data from paper [16] . Authors of this work applied modulated Ion Cyclotron 
Resonance Heating with 50/50 duty cycle, a modulation amplitude of approximately 80% 
and different frequencies of modulation /. As the first step, we will use experimental 
data from [16] for plasma ion thermal diffusivities ~ (d\nA m /dR)~ 2 , n a (f) ~ 
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(da/dR)~ 2 for R « 3.4 m determined from modulated temperature amplitude A m and 
phase shift a, respectively, where R is a linear coordinate in plasma along tokamak major 
radius. Authors of cited work found that < K a . We will show that (a) the difference 
between ha and K a is a result of finite speed of the perturbations propagation in plasma, 
(b) it can be fully explained by the local time concept, (c) only = k is the real 
incremental thermal diffusivity. 

In JET tokamak plasma vertical minor radius is 2.1 m, while plasma horizontal minor 
radius is 1.25 m only. Therefore along plasma major radius temperature field can be 
considered approximately planar and one dimensional as it is usually assumed [17] . So 
we can apply periodical solution of modified according to the local time concept one- 
dimensional Fourier equation from our previous works [2], [3] 



A m {R) = A expl-Z (R - R )} 



(87) 



a(R) = - 



£ (R-Ro)+u(t p 



(88) 



\/3 fuJ 



(89) 



where k, A m , Aq, Ro, u> = 2nf denote incremental thermal diffusivity, modulated tem- 
perature amplitude, reference temperature amplitude, reference R coordinate and cyclic 

frequency, respectively; n = 0, 1,2,3 Formulae (87)-(88) prove that ka = K , and only 

ka = k is real incremental thermal diffusivity which is present in modified and classical 



Fourier equations. Formula (88) predicts that in periodical oscillating temperature fields 



perturbation traveltime rp includes an arbitrary constant n/f, while speed of thermal 
perturbations Cx = \dTpfdR\~ 1 does not depend on it. 



Differentiation of formula (88) brings 



da I ' co \ 
dR = V + 



(90) 



where cr = (drp/dR)^ 1 is characteristic perturbation speed, and £ Q = y/3/2 ^Ju/n a . 
When in plasma thermal perturbations propagate outward, there is Cr = Cx > 0. In 
the case of inward propagating perturbations there is cr = — cx < 0. Therefore only at 
infinite speed of the perturbations propagation In reality, when always 

speed of the perturbations propagation is finite there is ka = k ^ K a . 



From formula ( 90 ) follows 



cr 
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(91) 
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Figure 6: Experimental values of characteristic perturbation speed cr of ion plasma com- 
ponent in JET tokamak vs. u according to data from Ref. 



We see that if n a > ka, as measured by authors of cited work [16], then characteristic 
propagation speed is negative cr < (Fig. 6), and the perturbations propagate inward. 
Fig. 6, built according experimental data from cited paper, shows that speed of the 
perturbations propagation increases when frequency of temperature oscillations grows. 
The data also confirms theoretical prediction that in magnetically confined fusion plasma, 
speed of the perturbations propagation perpendicular to the field is small. Thus, we 
can conclude that analysed experimental data from Ref. [16] support the theory of the 
perturbations propagation and the local time concept. 

The cited paper contains also results of simultaneous measurements of electron and 
ion components temperature amplitude and phase, though authors give experimental 
errors only for ion component. We use these experimental data for determination of si- 
multaneous characteristic propagation speeds c in ion and electron plasma components 
for / = 8 Hz, 3.2 m < R < 3.6 m. From experimental amplitude distributions A(R) 



were calculated by least squares and formula (87) incremental thermal diffusivity, ref- 
erence temperature amplitudes and reference R coordinates n e = 1.315 m 2 /s, K{ = 
0.461 m 2 /s, A 0e = 200.6 eV, A 0i = 165.6 eV, R 0e = 3.07 m, R 0l = 3.18 m for elec- 
tron and ion plasma components, respectively. Then formula (88) with fixed n = was 



applied to experimental phase distributions a(R) and were calculated experimental per- 
turbation traveltime distributions Tp(R) and, by numerical differentiation, experimental 
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Figure 7: Simultaneous space distributions of experimental characteristic perturbation 
speed cr in ion (triangles) and electron (diamonds) plasma components of JET tokamak 
according to data from Ref. [T6] . 



characteristic propagation speed distributions cr(R) = (drp / dR)' 1 for ion and electron 
plasma components (Fig. 7). We see that (a) speed of thermal perturbations propaga- 
tion in JET tokamak perpendicular to magnetic field is small and depends on R, most 



probably because of varying magnetic field B (R) as it follows from formula (72); (b) for 



electron and ion components experimental characteristic propagation speeds cr(R) have 
the same negative sign; (c) determined values of experimental characteristic propagation 
speed cr(R) for ion component are close to experimental characteristic propagation speed 
for electron component. Therefore we can conclude that results of this check support the 
local time concept and do not contradict our theoretical suggestion, that speed of the 
perturbations propagation and the thermal perturbations traveltime space distributions 
in the electron and ion plasma components are the same. 



8 Conclusions 

1. Developed general statistical theory of the perturbations propagation and the local 
time concept are valid for all non-equilibrium systems in gases, liquids, solids and 
plasma with any macroscopic velocities. 
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2. Finite speed of the perturbations propagation is predicted for all transport pro- 
cesses. It depends on the system properties, the process type and macroscopic ve- 
locity component normal to the boundary separating perturbed and non-perturbed 
subregions. 

3. Propagation of perturbations in transport processes is described by primary waves 
solutions of derived eikonal type equation for the perturbations traveltime. 

4. Developed in the work local time concept allows to formulate modified kinetic, 
conservation and governing equations describing real transport processes with finite 
speed of the perturbations propagation. 

5. At the same temperature, speed of the perturbations propagation in ideal gas is 
the highest for energy transport (thermal conductivity) at constant pressure, it is 
the lowest for momentum transport in isothermal viscous flow with zero normal 
macroscopic velocity component, and it has intermediate value for mass transfer in 
isothermal self-diffusion. 

6. For metals, speed of thermal perturbations propagation is the same in Fermi elec- 
tron gas and lattice. At low, ambient and moderate high temperature, it is, approx- 
imately, half of Fermi velocity. In dielectric crystals, speed of thermal perturbations 
propagation is equal to speed of sound in the crystal. 

7. In plasma, speed of thermal perturbations propagation is the same in electron and 
ion components. In a case of strong enough external magnetic field, speed of thermal 
perturbations propagation perpendicular to the field is small. 

8. Published experimental data for femtosecond laser heating of thin gold films is in 
satisfactory agreement with modified coupled Fourier governing equations, resulting 
from the local time concept, and calculated theoretical speed of thermal perturba- 
tions propagation in metals. 

9. Published results of power modulation experiments in JET tokamak are in sat- 
isfactory agreement with modified Fourier governing equation, obtained from the 
local time concept, theoretical prediction of small thermal perturbations propaga- 
tion speed perpendicular external strong magnetic field, and support our theoretical 
suggestion that speed of the perturbations propagation and the thermal perturba- 
tions traveltime space distributions in the electron and ion plasma components are 
the same. 
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